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Abstract. Let -. R'* > C be a compactly supported function which satis- 
fies a refinement equation of the form 

</3(a;) = ^ Cfc(/3(Ax - fc), Cfc e C, 
fcsA 

where V C K'' is a lattice, A is a finite subset of F, and A is a dilation matrix. 
We prove, under the hypothesis of linear independence of the F-translates of 
ip, that there exists a correspondence between the vectors of the Jordan basis 
of a finite submatrix of L = [c^i_j]i .jgp and a finite dimensional subspace 
"H in the shift invariant space generated by ip. We provide a basis of "H. and 
show that its elements satisfy a property of homogeneity associated to the 
eigenvalues of L. If the function has accuracy ft, this basis can be chosen 
to contain a basis for all the multivariate polynomials of degree less than k. 
These latter functions are associated to eigenvalues that are powers of the 
eigenvalues of A^^ . Further we show that the dimension of W coincides with 
the local dimension of and hence, every function in the shift invariant space 
generated by kp can be written locally as a linear combination of translates of 
the homogeneous functions. 



1. Introduction 

Let r be a lattice, i.e. F is the image of Z"* under any nonsingular linear trans- 
formation and ^ be a dilation matrix associated to F, (i.e. ^(F) C F and all 
eigenvalues of A satisfy |A| > 1). We will say that a compactly supported function 
if : R — > C is refinable with respect to A and F, if it satisfies the dilation equation 

^{x) = J2('k(p{Ax~k), xeW^, (1.1) 

feSA 

for some finite subset A C F, and coefficients e C. 

The Shift Invariant Space (SIS) generated by ip is the space 

Si^) - |/:M'^^C:/(x) = ^yfe(^(x-f A:), G C, A: G F I . 
I fcer J 

Note that since (p is compactly supported, the right hand side of the previous 
equation is well defined. Even though the SIS is an infinite dimensional space, the 
fact that its generator is compactly supported yields "locally" a finite number of 
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generators. More precisely, let C M be a fundamental domain for the lattice F, 
and 

Ei'p) = U/e -.E—^C: f/E{x) = fix) Vx&EJe S{^)}. (1.2) 

The space E{lp) is finite dimensional. A canonical set of generators is the set 

{ip{x - k)/E ■■ E — >C; k such that |(Supp(v5) + k) n E\ > 0}. 

The algebraic dimension of the vector space E{ip) will be called the local dimension 
of S{ip) and a basis of E{(p) a local basis for S{ip). 

If if satisfies a refinement equation one may not explicitly know the func- 

tion, but from properties of the coefficients of 11. 1() one often can deduce (proper- 
ties) of these finite generators. One question is if it is possible to obtain a different 
set of generators with specific properties. In particular it is important to know if 
some (or any) of these generators can be chosen to be polynomials. The accuracy 
of if is the maximum integer n, such that all polynomials of degree less or equal to 
K — 1 are contained in S{ip). Hence, if ip has accuracy k, one can choose a local 
basis containing = X)s=o '^'^ linearly independent polynomials, where dg is the 
number of linearly independent monomials of degree s. 

In the one dimensional case, with dilation 2, the generator ip satisfies 1^9(2;) = 
Sfc=o Cfc¥'(2a; — k). In this case — n and the accuracy of ip is related to spectral 
properties of a finite matrix T. 

Precisely, under the hypothesis of linear independence of the integer translates of 
the generator, ip has accuracy k, if and only if {1, 2^^, 2^'^'*^^^^} are eigenvalues 
of the (iV -|- 1) X (iV -I- 1) matrix T defined by T = {c2i-j}ij=o,...,N, (the scale 
matrix) and there exist polynomials po, of degree(pi) = i such that each 

of the the vectors Vi = {pi(k)}k=o,...,N is a left eigenvector of T corresponding to 
the eigenvalue 2'' f |Dau88j . |CHM98j V Here the fact that the powers of 1/2 are 
eigenvalues of T is related to the dilation factor 2 in the dilation equation. 

So, if ip has accuracy k, we know k linearly independent functions in E(ip). This 
set of functions can then be extended to a local basis of S{ip). A natural question 
is, in which way can this completion be done? If some eigenvalues are associated 
to nice local basis functions, (i.e. local polynomials) would it be possible to extend 
this set of nice functions to a basis of E{(p) using the remaining eigenvalues? What 
properties will these new functions have? 

Blu and Unser |BU02| in the study of radial basis functions, and later Zhou |Zho02| 
gave the first clue for the answer to these questions. They showed that associated 
to an arbitrary eigenvalue A of the matrix T there is a function in S{ip) that satis- 
fies that h{2x) = jh{x). In particular, the monomial x'^, satisfies this property, for 
A = 2"*^. So the set of functions associated in this way to all the eigenvalues of T 
is a linearly independent set. Then in the case that the matrix T is diagonalizable 
it is possible to complete the local basis for S{ip). 

In |CHM05j the problem was completely solved. They showed that to each vector 
from a basis that gives the Jordan form of the matrix T it is possible to associate 
a function h in S{ip) that they called (A, r)-homogeneous and satisfies that 

{D2 ~ Xlfh = 0. 

Here D2/ is the dilation operator defined by 1)2/(2;) = /(§), and A is an eigenvalue 
of T. In particular, the functions associated to eigenvectors are (A, 1) -homogeneous, 
and correspond to the one's obtained before. These (A, r)-homogeneous functions 
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in S{ip) are linearly independent and provide a local basis. This local basis contain 
all the monomials x'^ within the accuracy. The generator (p can be completely 
obtained from this local basis. 

The goal of this paper is to carry on this study to M'^, with a general dilation 
matrix and an arbitrary full rank lattice. 

When moving to higher dimensions, the situation turns much more complicated. 
In analogy to the one dimensional case, we will consider functions that satisfy the 
relation h{Ax) — jh{x) (see section 0)). Here h : ^ C, A is a d x d invertible 
matrix and A G C. More in general, we will consider functions satisfying that 
{Da — Xiy h = 0. To avoid any ambiguity, we will say that functions satisfying this 
equation are in the class Ti.{A, A, r), in place to use the word homogeneous, since we 
will also be dealing with polynomials that are homogeneous in the standard way, 
(i.e. a polynomial p of degree s is homogeneous, if p(ax) = a''p{x), x £ M'* Va e R). 
Note however, that with this definition, for d = 2 the monomial h{xi,X2) — xiX2 
will be in H{A, A, 1) only if A is diagonal and A = -^^^a^- 

When trying to extend the notion of accuracy from one to higher dimensions, it 
became apparent that the fact that ip has accuracy k is not immediately related to 
spectral properties of a finite submatrix of L = [cAi-j]i,jer- The relation is much 
more subtle and involved (see .CHMQ^ . .CHMOOj . |(^HM00j V 

In spite of this, however, in this paper we are able to obtain an analogous result to 
the 1-dimensional case. Again we are able to show that a local basis of S{(p) can be 
obtained using solely functions from Ti.{A, A, r), where A is an eigenvalue for a finite 
submatrix T of L. The result is very pleasing, in the sense that the 1-dimensional 
results are completely recovered and one obtains a different way of writing the 
functions of S{ip), namely, each function in S{ip) can be written locally as a linear 
combination of the translates of functions in H{A, X,r) (cf. equation (|4.10|l 'l. 

In particular, if ip has accuracy k, then we will find linearly independent 
polynomials that are in the class 7i(^, A, r) for some eigenvalue A of T, and some 
r e N. 

The difficulty here is to find the appropriate matrix T. Since we are in M.'^, the 
indexes vary along a d-dimensional lattice, so to write L as a matrix, one has to 
order the points. Which order is not important, as long as it is always the same. 
In the one dimensional case, it was straightforward to look at a submatrix of L 
that was intimately related to the support of tp. In the higher dimensional setting, 
it may be a difficult problem to determine the support exactly. This is one of the 
problems one has to overcome to solve the question raised here. 

The paper is organized as follows: In section|21we briefly review some geometric 
properties of the support of a refinable function tp related to the dilation A and 
define the finite matrix T whose spectral properties will be fundamental for our 
analysis of the class 'H{A, A, r). In section|21 we relate the spectral properties of the 
infinite matrix L to those of T, and in section0]we define the class Tl{A, A,r) and 
show how one associates one of these functions to each vector of the Jordan basis 
of T. We further prove that these functions are a basis for the space E{(p) H1.2|l . 
Finally, in section O if (/j has accuracy k, using results from CIIM03 , we show 
that the space of all functions in the class TL{A, A, r) in S{(p), contains ak linearly 
independent polynomials. 
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2. Attractors, Tiles and Admissible Sets 

Let r be a lattice and A a dilation matrix associated to T. Then A has integer 
determinant and the group r/A{T) has order |det(A)| (see for example |Woj97| ). 
Set 

TO = |det(A)|, 

and let D ~ {di, . . . , c?™} be a set of representatives of the group r/A{T) of order 
TO. We call D a full set of digits, or digit set. 
The cosets 

Ti = ^(r) -d, = {Ak -di-.ke T}, d^ G D, 

form a partition of F. Assume that 71, . . . , 7^ is a set of generators for F, that is, 
71 , . . . , 7d are linearly independent vectors in R*^ and 

r = {iai + ... + ia-fd-heZ}. 

We will call the set 

P = {xi-fi + . . . + Xdjd : < a;^ < 1} 
a fundamental domain for the group R'^/F. 

2.1. Attractors. For each A: e F, we define Wk : R'' — yMf^ by 

Wk{x) = A~''^{x + k). 

Since A is a dilation matrix, A^^ is contractive for some appropriate norm in R"*, 
so each Wk is a contractive mapping on for that norm. 
The space 

niM"^) = {K : K ^9 and K is compact}, 

is a complete metric space under the Hausdorff metric d defined by 

d{B, C) = inf {e > : S C and C C B J , 

where 

= e R'' : dist(a;, B) < e j . 

For each finite subset H cT, we define 

WHiB) = U WkiB) = A-\B + H) 
ken 

It can be shown that wh is a contractive map in Ti.{W^) (using that each Wk is a 
contractive mapping on R'^). Consequently, by the Contraction Mapping Theorem, 
there exists a unique nonempty compact set Kh C M.'^ such that 

WHiKH)^KH, I.e. Kh^A-\Kh + H). 

In fact, we can write 

CSO I 00 I 

Kh^Y. ^"'(^) = { E -hjeHy (2.1) 

The set Kh is called the attractor of the iterated function system generated by 

{wk}keH. EntHn. 
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2.2. Tiles. Given tiie digit set D — {di, . . . , dm}, we consider the attractor 

oo I oo I 

Q = Kn = Y^ A-^{D) = I -ejeDS (2.2) 

j=i [j=i J 

of the iterated system generated by {wd]deD- We have that, for 7 G F 

00 00 
KD+-, = ^"'(^ + = ^"'^ + " ^y^^' = Kd + {A- (2.3) 

Therefore we can assume without loss of generality that G -D, and hence, by 
equation H2.2|l . we have G Q. 

The set Q satisfies the following properties (see jBan91| and |GM92| '): 

a) [jkerQ + k=^''- 

b) 0,<9 = Q°, and \dQ\ = 0. 

c) IQ n (Q + fc)| = for every k eT - {0} if and only if \Q\ = \P\, where P 
is a fundamental domain for W^/T. In this case Q Cl {Q + k) C dQ for all 
k e r - {0}. 

A longstanding problem was the question of whether for each dilation matrix 
A there exists a full set of digits D such that the corresponding attractor Q is a 
tile. A counterexample was found recently (See |Pot97j and also |LW99. .) We will 
assume in this paper that Q is a tile, or a fundamental domain for F, that is, the 
F- translates {Q + fc}fcgr cover with overlaps of measure zero (hence \Q\ = \P\). 
Then the local dimension of S{lp), will be the dimension of Q{(p). See (|1.2|l . 

2.3. Admissible sets. Let _ff be a fixed finite subset of F 
Definition 2.1. We say that a set £7 C F is H- admissible if 

A-\n + H)nr cn, (2.4) 

which is equivalent to say that wh{^) H F C fi. 

Remark, li H C H and H. is H -admissible, then fl is i?- admissible. 

We immediately have the following Proposition: 

Proposition 2.1. If fin is defined as fin — Kh H F, then fin is an H -admissible 
set. 

Proof. Since ^lH C Kh, we have 

WH{nH) n F C Wh{Kh) r\T^VLH, 

which shows the desired property. □ 

Let ^(F) be the space of all sequences defined in F, and let L be the infinite matrix 
associated to the refinement equation Lij — CAi-j, whenever Ai — j G A, and 

Lij = in all other cases. 

In this paper we will mainly be interested in A-admissible sets. The reason for 
that is that if $1 C F is A-admissible, then the space £{fl) — {Y G ^(F) : yk — 0, k ^ 
n} is right invariant under L. 

We will need to "extend" finite vectors to infinite ones with certain prescribed 
properties, and such that they coincide with the finite one if restricted to a finite 
subset of the lattice. Therefore, the following Proposition found in |CHM04j . will 
be very useful. 



6 



C.A.CABRELLI, S.B.HEINEKEN, AND U.M.MOLTER 



Proposition 2.2 (CHM04). For each finite H <zV, there exists a strictly increasing 
sequence {rin}n>o of H- admissible sets whose union is T, such that fJo = and 

ii'H(f^„+i)nrcr!„ (2.5) 

for all n > 0. 

Proof. Let || ■ || be any norm in M'* such that ||^^^|| < 1 and fix e > 0, such that 
H C B{e), where B{e) = {x G M'' : ||a;|| < e}, the closed ball with radius e centered 
at the origin. Now set 

Choose i5 > 5o in such a way that C B{S) and set Fq = B{S). Since S > So, we 
have ||A"i||((5 + e) < S. Hence, 

wh{Fo) = A-\B{S) +H)c A-\B{5 + e)) c B{\\A-^\\{5 + e)) c B{5) = Fq. 

We define recursively F^+i = wniFj) for j > 0. It is easy to see, by induction, that 
Fj^i C Fj for every j. Since Fq is compact, the Contraction Mapping Theorem 
tells us that P| Fj = K^. It follows that Fj OT = fiji for every j large enough, and 
consequently {Fj n F} is a finite collection of sets. Let 

f^ff = f^o C f7i C ■ • • C f^AT = Fo n F 

be the distinct elements of this collection and fix < n < A^. Since there exists a 
j e N such that 

f}„ = F, nFcFj_inF = r!„+i, 

we have 

(f^„+i) n F c wh{Fj-i) n f = n f = f7„. (2.7) 

So inclusion l|2.5|) holds for rt = 0, 1, . . . , iV - 1. 

Now we set 6n = S and define recursively, Sn+i = — e for n > N. The 

sequence of numbers Sn < Sn+i < ■ ■ ■ is increasing. Define il„ = B{6n) H F for 
n > N. If f2„+i = f2„, we skip that one and continue until Qn+k ^ ^n- In this way 
we obtain a strictly increasing sequence of sets {r2fc}fc>Ar. Combining with the sets 
ilo, . . . , r^AT constructed previously, we have a strictly increasing sequence {r2n}„>o- 
The inclusion 

WH(f^„+i)nFcr!„ (2.8) 

holds for every n G No, since for n > iV, again there exist a j e N such that 
0„ = B{8j) n F C n F = r!„+i, 

and then 

WH{^n+x) - A-^rin+l C B( 1 1 A" ^ || ( + i + e) ) B{8j). (2.9) 

We already showed that fio = i^ff is 7?-admissible. Since fin C iln+i, it follows from 
(j2.8(l that is if-admissible for every n e No, which completes the proof. □ 

Corollary 2.1. If H C. H' d T, then there exists uq > 1 and a strictly increasing 
sequence {^ln}n>o of H- admissible sets whose union is T, such that 

i7o = ^H, — ^H', o,nd wni^n+i) n F C fin for all n > 0. 
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Proof. First construct the sequence {i^'n} associated to H' using the previous propo- 
sition. Note that by Remark 12.31 the sets fi^ are also i/-admissible. 

Using the notation of the previous proof, let jo be such that Fj^ fiF = fin'- Now 
consider the sequence {Gj}j>jo, where Gj„ = Fj^ and Gj+i = wniGj), for j > jo- 
Since 

G,„+i = WHiG,„) = A-\F,„ +H)C A-\F,„ + H') = wh'{F,„) C F,, = G,„, 

then Gj+i C Gj and therefore r]j>jgGj — Kh and hence, {Gj nF} is again a finite 
collection of sets of say no + 1 elements. Consider now the distinct elements 

riff = r!o c r!i c ... c = n f = ri^,, 

and let 

This new sequence satisfies all the desired properties. □ 
For a more complete treatment of admissible sets see |CHM04j and also Jia |Jia98 | . 

3. Spectral Properties of L 

Let us now return to the refinement equation (p{x) = X^feeA '^kfiAx — k). 

If we consider the infinite column vector 

= {(^(x + fc)}fegr, (3.1) 

this equation becomes 

^{x)=L^{Ax). (3.2) 

It can be shown (see jCHMOOj ) that the set K/^, which is the particular case taking 
_ff = A in H2.1(l , satisfies that if is a compactly supported solution of the refinement 
equation then Supp((^) C i^A. Also, by Proposition l2.1l the set 17a — ^aHF 

is A-admissible. However, it is not necessarily true that Suppc/j C UagOaQ + 

We will therefore consider the bigger set SI' — K\r n F, where A' — A — D D A. 
In KjhM04| it was shown that the translations of Q using all elements of fi' cover 
the support of the compactly solution to Hl.l|l . Moreover, SI' is A' admissible, and 
hence also A-admissible. As noted earlier, the A- admissibility of il' guarantees that 
the space £(S1') = {Y S ^(F) : yk = 0, k ^ fl'} is right invariant under L. 

Let now {S7„}„>o be a sequence of subsets of F that satisfies: 

• f2o — ^A 

• For i > 0, fli C S7i_|-i, and Ui^li = F 

• For z > S7i are A-admissible and WA{^i+i) n F C fi;. 

• n^o = ^' 

• For i > uq, fli are A'-admissible and WA'(S7i+i) HFC fj^. 
These sets exist by Proposition 12 . 21 and its Corollary 12. II 

We denote by {T'„}„>o the finite submatrices of L 

Tn = [cAi~j]i,jen^- (3.3) 

Since S7„ C S^n+i, if the order in F is appropriately chosen, actually r„ is a subma- 
trix of T„+i, for each n. 

Let Y = {yk}k£T £ ^(r) be an infinite row vector, and P„ : ^(F) — > (j^ixo,,^ ^ ^ 
be the restriction mappings defined by 

PnY^{Vk}km^- (3.4) 
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We consider L — XI : £{T) — > ^(r) the left-multiplication operator who maps 
Y — > y{L — XI) (where / is the identity operator acting on ^(F)). By abuse of 
notation, / will be any identity operator, no matter on which space it is acting on. 

Note. In what follows we will use powers of the matrix {L — XI). Note that these 
powers are point-wise well defined, since the rows of the matrix L have a finite 
number of non-zero elements. 

The next proposition shows the relation between the spectrum of L and r„: 

Proposition 3.1. Consider A £ C, r G N and n > 0. 

(1) Let Y € £(r). We have 

Y e Ker(L - XI)'' implies P^Y e Ker(T„ - XI)\ (3.5) 
Conversely, 

(2) If V £ Ker(T„ — XI)'' and A 7^ 0, then we can extend v to an infinite row 
vector Yy (i.e. Yy £ £{T) and PnYy ^v), so that Yy G Ker(L - Xiy . 

(3) If X ^ 0,Y and Y £ Ker(i - XI)'', then P^Y ^ 0. In particular the 
extension in ^ of v to Yy is unique. 

Proof. 

(1) First note that j G fin and Ai — j G A implies i G fin- For, in this 
case, Ai G ri„ + A and since iln is a A-admissible set it follows that i G 
A"i(ri„ + A) n F C 0„. Hence 

j G nn, I ^ =^[L- XI],j = 0. (3.6) 

Moreover, we will show by induction on r that, 

if j G rin and i ^ rj„, then [(i - XI)'']^. = 0. (3.7) 

(a) The case r = 1 is simply (|3.6I) . since we assume that = if fc ^ A. 

(b) Suppose now that (|3.7|l holds for some fixed r > 1. Using (a), for 
j G r2„ we have 

ker 

= E [iL-Xm^.iL-XI]^,. 
fceo„ 

Now, if i ^ fin, the inductive hypothesis yields that the last sum is 
zero. 

Therefore, the statement is true for all r G N. 

To prove the first part of the Proposition, let Y G ^(F) and Y G Ker(_L — 
XI)'-. 

Applying the preceding equality, we obtain for each j G On 

[{PnY){Tn-xm^^j2y^i(^-^^yh 



Ey»[(i-A/)1,,. 
[r(x-A/)n, = o. 



This completes the proof of (1). 
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(2) Assume that v G C^''^^",A / and v £ Kcr(T„ - \iy . Wc want to 
construct a vector Yv G £{T) such that P„K„ — v and Yv G Ker(L — Xiy . 
We now prove by induction on r that, 

for J e n^+u * ^ f^fc, fc > 0, [(L - Xiy] = r f ' ' ^ ^. (3.8) 

I (~^) lo'" * = J- 

(a) Case r = 1. If z were such that Ai-j £ A, then i e A^^(rife+i+A)nr C 
rifc. Hence, ii i ^ ilk then Ai — j A and so = 0, and therefore 
[L - XI],j = for i ^ j, and [(i - XI)]jj = -A. 

(b) Assume that (|3.8(l holds for r > 1. Then for j e ilfc+i and i ^ flk, k > 
0, we have 

[(L - A/)'-+i],^. = 5: [(L - A/)'-],jL - A/],, 

= ^ [(L - A/)'-],jL - A/],„ 

since by the inductive hypothesis [{L — Xiy]^(, = if ^ G fi/c C ^k+i 
and i ^ flk- It foUows using the case r = 1 that 

which proves (|3.8() . for every r G N. 
Define now ?/j = for j G and define recursively, for j ^ ri„. 



^- ^ y,[(L-A/)n„, J Gr!fc+i\r!fc,A:>n. (3.10) 



The vector y„ = {yj}j^r is an extension of v. To see that Yy G Ker(_L— A/)'', 
since (y(L - Xiy)j = J^ier vAi.^ - ^^YVoi we have: 
• If j G Vln, then by (j?^ 



^y.[(L-A/)'-],^.= 5^y4(r„-A/)n,^. 

= ^«,[(T„-A/)n,^. =0. (3.11) 

iGf2„ 

• If j ilji, then there exists k G No,/c > n such that j G ^k+i \ ^k- 
Therefore, 

- A/)'-],^. = 5^ - xn% + E - 

= J2y^[iL-xIy]^^+y,i-xy) 

= 0. (by ^J^) 

(3) For the last part of the Proposition, assume that A 7^ 0,y ^ 0, and Y G 
Ker(L - Xiy. To show that P„F ^ 0, take fco G T such that yk„ ^ 0. If 
ko efln, there is nothing to prove. Otherwise, let 

to = min{A: G N : fco G Qk}- (3.12) 
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Since Y{L - Xiy = 0, 

So, there exist fci G fi^, < fc < io, such that yk^ 0. If fci G fi„, we can 
stop here. If not, we repeat the procedure until kj G fin- 

□ 

Remark. 

• Since the previous Proposition is true for any set of the sequence fin, in 
fact the smallest matrix Tq already has all the spectral information of L. 

• The extension of the vectors of Ker(ro — A/)'' to vectors of Ker(L — Xiy 
will produce intermediate vectors of Ker(T„ — Xiy , by the construction of 
the sets f2„ produced in Corollarv l2.1l 

For the special case A = 0, under some mild assumptions, we have an additional 
property. We say that the F translates {ip{- — k)}ker are linearly independent, if 
for any sequence {ak}ker in ^(r), 



akip{- — k) — implies ak — 0. 



ker 

Lemma 3.1. // {Lp{- — k)}ker o,re linearly independent, then the operator L : 
^(F) — * t{V), r I — >YL is one to one. 

Proof. Let YL = Q. Then 

Y^{x) = YL^{Ax) = 0. (3.13) 

Since {ip{-—k)}ker are linearly independent, Y<^ — imphes F = 0, so Ker(L) = {0}. 

□ 

4. The class 'H{A,X,r) 
Assume Y G Ker(L — Xiy . If we define h{x) ~ we have 



Q = Y{L- Xiy^x) = Y {-Xy-^L^j ^{x) 

= ^ (E (I) ("A)'-'=$(A-^x) j , 



fc=0 

This leads to the following definition: 

Definition 4.1. A function h is in the class T-L{A, A, r), if it satisfies 



k=0 



E ( I ) {~>^fHA^^''x) ^ for every a; G M''. (4.1) 
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If we define the operator Da by DA[f )[x) ~ f{A ^x), then h is in H{A, A, r) if 
and only if 

{Da - Xiyh = 0. 

A function in 7i(^, A, r) will also be said to be of class 7i(^, A, r). 
Note that if /i G 'H{A, A, r), then h € 7i(j4, A, s) for every s > r. 

Proposition 4.1. Let C M'' 6e a bounded set such that E V° and V C AV. 

Set C — AV \ V and a G Z. Let h be a function of class Ti.{A, A, r) . Then the values 
ofh in R'' \ {0} can be determined from its values in any set of the type: 

c IJ A'^c. 

k—a 

Furthermore, if X I then h(0) — 0. 



Proof. Since h G Ti-iA, A, r) we get that 

r 
k 



k=l 



m = -Y.[ I )i-^)''HA'x) and (4.2) 



Hx) = -T.( I )i-^)-'h(A-'x). (4.3) 



k=l 

On the other side, it has been proved in |ACM04j that the set C satisfies: 

a) U,eZ^'C = R''\{0} 

b) The sets {A^C}^^'^ are pairwise disjoint. 

So, from H4.2() we deduce that the values of h in A°~^C can be obtained from the 
values in C, and analogously, from (|4.3|l the values of h in A°'~^^~^^C can be obtained 
from the values in C. 

Then we proceed inductively to obtain all the values in R'* \ {0}. Finally, it is 
immediate from the definition, that h{0) = when A 7^ 1. □ 

Proposition 4.2. Suppose {(p{- — k)}ker o-te linearly independent. Let /i, . ■ . , /; G 

S{(p), fi = where V G £{T) . Then /i, . . . , // are linearly independent func- 

tions if and only if .... .Y'^ are linearly independent in £(r). 

Proof. Since 

; I / ' \ 

J2 = E "'(^'*) = E (4-4) 

i=l 1=1 \i=l / 

and the translates of (p along the lattice F are linearly independent, we conclude 
that J2\=i^ifi = if and only if {J2\^iCeiY^) = 0, which leads to the desired 
result. □ 

Remark. Let E : S{(p) — > £(F) be the function that associates to each element of 
S{(p), its coordinates in {(p{x+k)}. Proposition l4.2l shows that E is an isomorphism. 

Proposition 4.3. Assume that {ip{- — k)}ker Cire linearly independent. 

(1) If h & Siip),h = F$ and h G HiA,X,r), then Y G Ker(L - Xiy and 
P„r GKer(r„-A/)''. 

Conversely 
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(2) Assume that A ^ 0,f G Ker(T„ — Xiy and that Yy is the unique extension 
of V such that Yy G Ker(L — Xiy (see Provosition \S. Then the function 
/i = 1;,$ belongs to H{A, A, r). 

Proof, li h = Y^ is of class 7i(A, A, r), then we have 

r 
fc=0 

= E (l) (-A)'=rL'-'=$(.T) 

= y(L- A/)''$(a;). 

Since the F translates of (/? are linearly independent, it follows that Y{L — Xiy = 0, 
and consequently, by Proposition Kill PnY G Kcr(T'„ — Xiy. 

To prove the second part, note that if v = 0, the statement is trivially true. If 
A 7^ and v G Ker(T!„ — A/)'' and v ^ 0, then, by Proposition l3.il we can extend v to 
a vector Yy G Ker(L — Xiy , and so the function h — Yy^ is of class 7Y(A, A, r). □ 

4.1. Jordan decomposition of T„. Let to„ = Consider the set A„ of 

eigenvalues of T„ and the associated Jordan basis i3„ = {vi , . . . , } of C"" . For 
each Vi G y8„ we have that Vi G Ker(T„ — Xiy and 1;^ ^ Ker(T„ — A/)'^~^ for some 
A G A„, and for some fc > 1. So to each Vi there corresponds a unique pair (A, k). 
Note however that to two different ViS in the basis there could correspond the same 
pair. For each vector of Bn, set Vi = Vi{X, k). If A ^ 0, by Proposition 14.31 we can 
associate to each Vi{X, k) a function hy.(^\ i^-) in 7i(A, A, k) f^S{^p). Since the vectors 
Vi(X,k) are linearly independent, its extensions {i^;} are linearly independent in 
^(F), so the functions {hy.(^\^k)}vieBri,\^o are linearly independent. 

If hi,..., hi arc of class Ti.{A,X,ki), for some ki,i = 1,...,Z, then a linear 
combination of them is of class Ti.{A, X, k), with k = maxi<i<;(fci), for if h = 
J2\=o(^ihi{X,h), then 

/ I 

{Da - Xlf h^{DA- Xlf a,h, = ^ a» {Da ~ XI f h, = 0, 

1=0 i=0 

and consequently, h G H{A, A, k). 
Let 

XtJx)= n i^-^Y' (4-5) 
AeA„ 

be the characteristic polynomial of T„, and set 

H\{(p) = {he S{ip) : h G H{A, A, fc), for some fc > 1 }, A G A„. 

Then, if A 7^ 0, using Proposition ^31 <i^^{'H\{'f j) — r\ and a basis for H\{ip) are 
the functions of class H{A,X,k) corresponding to the vectors Vi G Bn, such that 
Vi = Vi{X, k), for some A; > 1. So, if we denote 

n= nx{^)cs{v), 

then, dim(?i) = mn — ro, where tq is the dimension of the subspace generated by the 
vectors of the Jordan basis associated to A = 0. (Note that m„ = X^asA^ a^^o'^a-) 



n {~xyh{A''-^x) 



REFINABLE SHIFT INVARIANT SPACES IN 



13 



In order to be able to include the case A = in our analysis, we need to consider 
the case in which Supp(/3 C Ucjgo„Q + ^- This will guarantee, that (except for a 
possible set of measure zero), (p{x + /c) = if /c ^ 

4.2. The case in which ri„ contains Supp((^). If we recall the choice of the 
sequence {^n} at the beginning of section|31 it is clear, that for n > uq, we have that 
Supptp C Uujen„Q + ^, and hence, the local dimension of S{ip) is dimspan{(^(a: + 

In that case for x G Q°, X ^ 0, /i„i(A,/c)(a^) = + since 

ip{x + i) = if j ^ f7„. 

Moreover for A = we have the following Proposition: 

Proposition 4.4. Let n > uq, and let rg be the power of x in xt„ (c-f- (|4.5|l ). 

Consider Vi — Vi{0,r), with r < rg. Define h{x) — J2ken b^ilk^ix + k). Then 
h = a.e. on Q. 

We postpone the proof to remark that with this Proposition, if rt > uq, and Bn 
is (as before) the matrix whose rows are the vectors of the Jordan basis for T„, then 



(x) 



BnPn^ix) a.e. X € Q. 



(4.6) 



Hence, since the matrix Bn is invertible the local dimension of S{ip) coincides with 
dimspan{/i^.(_\ ^.^(x), a: G Q,Vi G Bn}, which is equal to the dimension of Ti.. 

So the local dimension of S{ip) can be found by finding the Jordan form of any 
of the finite matrices T„ as long as n > tiq. 

Moreover, any function of the shift-invariant space S{ip) can be written as a 
linear combination of the lattice translates of the homogeneous functions. Namely, 
let / e S{ip), then 

/(a;) = ^ a^(p(x + 7) xeR'^,a^eC. (4.7) 
If we call g : M" — > C"" and h : M" — > C"'" the functions 



g{x) = Pn'^{x)xQix) and h(a;) 



h. 



■vi{\.k) 



(x) 



,fc)(a;) 



Xq{x), (4.8) 



and for 7 € F we denote by = (aij+^y, . . . ,ai„^+-y) the vector of length m„ 
whose indices are in ri„ + 7 (here = {ii, . . . , im„}), then H4.7|l becomes 

/(x) = ^S^g(x + 7), (4.9) 

and using 1)4.6(1 we obtain 

f{x) = p^h{x + 7) where f3^ = a^Bn^. (4.10) 

We will now prove Proposition 14.41 For this, let rg be the power of x in xt„ 
(c.f. (|4.5|l '). Choose m > n large enough such that 

nm^rtn,- {D + AD + --- + A'"-^D) . (4.11) 
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Define the matrices [{Tk)d]ij = CAi-j+d,hj G ^^fe, for any k eN, and d E D. It is 
sliown in |CHM04j . tliat ii x E Q, for any r > 1 there exists G Q, 7r € T, such 
that 

X = A^^'iur + 7r) with jj. = dr + Adr-i H h A^^^di (4.12) 

where di G D, 1 < i < r. Therefore, if fc > no and is as in (13. 4|) 

Pfe$(x) = {Tk)d, . . . {Tk)d,PkHA'-x -jr) xe Q. (4.13) 
For convenience, we will call Q — flm- 

Lemma 4.1. With the previous notation, if r < ra, then for fc G f2 and j G ^ng, 
we have 

[-^m]fe(j-7^) ~ [{Tm)di ■ ■ ■ {Tm)dr]kj ' 

where jr e D + AD -\ h A'-'^D. 

Remark. Note that the preceding equation does not state that both matrices are 
equal. 

Proof. We will prove the Lemma by induction on r. Let 7^ be as in (|4.12|) . 

• The case r = 1 is trivial by the definition of {Tm)di- 

• r — 1 r Observe first that by the choice of fi, 

[^m]«(i-7^) = [(^™)<ir]«(j-yl7^_i) i U G fi, j G ^uq- 

Now 

['^m]fc(i-7,) = X! t-^™ [-^™]«0-7,0 
tiGO 

= ^ [^m ] fc„ [(^"0<ir]«(j-A7^_i) 

= E [T:^\J(T,n)dXu-A,.^.). (4.14) 

tiGO„,-|-7r-i 

where the last equality follows from the A'-admissibility of rj„f, . But 

[{Tm)dAu(j-Ai^-i) — [(^'")'ir](«+7^_i)j J " ^ ^"o ~ lr-1, j G , 

and therefore, using induction and the A'-admissibility of , H4.14|l be- 
comes 

[^m]fc(j-7^) = E ["^m "^]fc„ [(^™)rf.-](,i+7^_i)j 

Mer2„f,-7r-i 



fen 

which completes the inductive step. 



□ 



We can now prove Proposition 14.41 
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Proof. Let n^ilmhe as before, and let x G (9\ (^QuULi A-'dQ + A-'D-\ h 

A^^D). Note that with this choice of x, A^x — 'jr ^ Q° , and therefore, if u ^ ^na, 
then ip{A^x — 7^ + u) = 0. 

Using this, together with the previous lemma, and equation (|4.13|l for k = m, 
we have 

[T;,P^^{A^x)], = ^(^"^^ + J') 

'J en 

= XI ■ ■ ■ - 7r + j) 

jen 

= f{x + k). 

Therefore, for x G Q \ {dQ U ULi A-^dQ + A-'Z? + • ■ • + A-^D) 
h{x) ^ 'Y [vi]k'.p{x + k) 

feGO,, 

fceo,, jsn 

= E( ENfera,J^(^'^^ + j) = o 

□ 



5. Accuracy and homogeneous polynomials 

In this section we will relate the previously obtained results, to the "accuracy" 
of a scaling function. We will use the notation of jCHM98| . 

Definition 5.1. The accuracy of <^ is the highest degree k such that all multivariate 
polynomials q with degree {q) < k are in S{ip). 

Let X — (xi, . . . ,Xd) G K''. With the standard multi-index notation we write 
, where a = (ai, . . . , a^) with each ai a nonnegative integer. 
Denote by \a\ = ai + ■ ■ ■ + ad- The number of multi-indices a of degree s is 

\s + d- 1' 
d~l 



ds = 

For each integer s > 0, we define the vector valued function X^^^ 



by 

X[,]{x) = [x"]i^l^s, xeR''. 

The ordering of the multi-indices a of degree s is not important as long as the same 
ordering is used throughout. 

We will now look at the behavior of X^^] (x) under the multiplication by an 
arbitrary d x d matrix Z with scalar entries Zij. If |q;| = s, then (Zx)" is not 
in general a monomial. Instead, it is a new polynomial of degree s, that is still 
homogeneous, but possibly involves all terms x'^ with = s. 
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Let = [<^^]|a|=s,|/3h 

defined by the equation 



be the ds x ds matrix whose scalar entries ^ are 



7 , ^a,f3^ 



|/3| = s »=1 

The matrices and their properties have been intensively studied in |CHM98j . 
|CHJV199j . In particular, if Id denotes the identity matrix in R'', we have that 

(Id)[s] = Id, 

and if Z and U arc two matrices, 

{ZU)[s] = -^[sj^s] and hence, if Z is invertible {Z^^)[g] = (Z^g-^)^^. 
Dilation of X^g] {x) by Z obeys the rule 

X[s]{Zx) = Z[g]X[g]{x), 
hence, if A is the dilation matrix corresponding to the refinement equation (jl.lfl . 

If Jg is the Jordan form of A^^j^], then there exists an invertible dg x ds matrix Qs 
such that QsA^^[s]Qj^ = Js- So we have that 

QsX^s^iA-'x) = {QsA-\s^Q-')QgX^s]ix). 

Denote by Qs{x) = QsX[s](x). Observe that Qs{x) = {q\{x) , . . . , Q'!," {x)^ is a 
column vector of polynomials of degree s that are homogeneous. By the previous 
equation, we have that 

Qs{A-^x) = JsQs{x). 
Let (3 be an eigenvalue of and B the Jordan block of order £ associated to 

/?, i.e., 

( 13 ... \ 



B 



1 (3 


\ 





(3 

1 p J 



We write Qb{x) for the vector that is the restriction of Qs{x) to the coordinates 
that correspond to the block B, i.e, if j, j + 1, . . . , j + ^ — 1 are the columns of B in 



Js then Qeix) = (Qi(x), . . . , Ql+'^-^x)'^ . Since 



QsiA-^x) = JsQs(x), 

we have 

Qb{,A~^x) = BQb{x). (5.1) 
This relation will enable us to show how, under the hypothesis of accuracy, we 
can relate the Jordan form of ^"""^[s] to the one of T„. This relation also gives a 
necessary condition for to have accuracy k. 

Proposition 5.1. Assume that Lp has accuracy k and that {ip{-'-k)}ker are linearly 
independent. Let s < k. If P is an eigenvalue of A^^i^s] and B is a Jordan block 
of A~^^s^ associated to (3 of order I, then Tn has a Jordan block associated to j3 of 
order f with i' > I. 
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Proof. Consider Qb[x) — {Q^g{x), . . . , Qg{x)). It follows from l|5.1|l that 

Q],{A-'x) = /3Qij(x) 
QliA-'x) = Qh{x) + PQl{x) 

(5-2) 

QUA-'^) = Q'B'\x)+f3Q'si^). 
Since Q^(a;) G for 1 < i < ^, we can write 

%(a;)=y^$(a:), 

for some infinite column vector Y^. From (|5.2|l we have for 2 < j < ^, 
Y''^{A-^x) = Y''-^^{x) + /3r'$(x), 

which implies 

Y'L^ix) ~ l3Y'<^{x) = Y'-^^{x). 
So, the linear independence of {(p{- — fc)}fegr yields 

Y\L - 131) ^Y'-\ (5.3) 

Since ^^(a;) G n{A,l3,e), by Proposition O we have that w = P„y^ G Ker(r„ - 
piy. Consider the vectors wi = P„r^ f2 = {PnY'^){Tn- (31), ...,ve = (P„r^)(T„- 
f3lY~^ . Let us show that wi, . . . ,ve are linearly independent: Assume that 

^ a,«, = 0. (5.4) 



Since 



J2 c^r^r (Tn " PI)'-' = ^ a,z;(r„ - /3/)'-i (r„ - pi)'-' 



\i=l 



e+i-2 

i=l 

aiviTn - PIY 



it follows from H5.4(l that 

aMTn - Plf-^ = 0. 
Since for every Y E i{T), r G N we have that 

(PnY){Tn - piy = Pn{Y{L - Plf), 

part 3 of Proposition 13.11 tells us that v{Tn — piY~' ^ 0. Hence Ui = 0. If 
we multiply each side of (|5.4|l by (T„ — pi)'~^ we see that 012 = 0. Analogously 
Qfa = . . . = = and therefore vi, . . . ,Vi are linearly independent. This implies 
that we have a Jordan block of T„ associated to P of order at least £. We can 
repeat this procedure for every Jordan block Bi, . . . , Bk of ^"""^[s] associated to P 
of respective orders h > h > ■ ■ ■ > h- Let, for 1 < j < k 

Q%(x)^yI^^x). 
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All we have to prove now is that 

p„r/S (p„y/0(T„ -/?/), (F„y/0(r„-/3/)'i-i, 

are linearly independent. Let 

alRaYl'+ al{P^Yl'){Tr, - m+ •••+ a'^ {PnYt){Tr, ~ [31^^ + 

alPaYt+ al{P^Yt){T^-(5I)+ ...+ (P„r,''=)(T„ - = 

(5.5) 

Let Bi, . . . ,Bt the Jordan blocks of order li. If we multiply each side of the previous 
equation by (r„ — piY'^^^, we obtain 

t 

Y^^liPnYhiTn-Piy^-'^O, i.e. 

1=1 

Pn(j2alYl^{L^piy^-'^^0. 

Since X;Li "j^/H-^-/?^)''"^ £ Kci{L-f3I), part 3 of PropositionlOimphes that 

t 

i=l 

So, since by (ESJ and Proposition Ol ^/^(L - /3/)'i~\ . . .,yI^{L - (31)'^^^ are 
linearly independent, it follows that a\ = . . . = a] = 0. Repeating a similar 
argument for every lj,2 < j < k we can see that every scalar of H5.5f) is equal to 
zero. This completes the proof. □ 

Let us now recall H5.1|l , and notice that 

Qb{A-^x) - PQb{x) = {B- I3I)Qb{x). 

Equivalently, if we recall the definition of Da of the previous section, DA{f){x) — 
f{A^^x), we have 

{Da - mQsix) = {B- (3I)Qb{x), 
where the product on the left side is understood coordinatewise. Moreover, for 

{B - pifQBix) = ( ■) {-f3)''-'B^QB{x) 

i=0 ^ ^ 
k 



i=0 



= {Da- pITQBix). 
In particular, since {B — (31) is nilpotent of order £, we have 
{Da - PIYQb{x) = {B- PIYQb{x) = 0. 
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Hence, all entries of Qb{x) belong to 'H{A,(3,£). We can repeat this argument for 
every Jordan block associated to f3 and every eigenvalue /3 of A^^^g]- It follows that 
each component of Qs{x) belongs to A, r) for some eigenvalue A of ^~^[s] , and 
some r S N. Since Qs is an invertible matrix and the monomials with \a\ = s 
are linearly independent, it follows that Q\{x)^ . . . , Q'^'{x) are linearly independent, 
and all homogeneous polynomials q{x) = q{xi, . . . , Xn) with deg{q) — s, are a linear 
combination of Ql{x),. . . , Qg'{x). 
We can now state the next theorem: 

Theorem 5.1. Assume that tp has accuracy k and that {ip(- — k)}ker are linearly 
independent. If q is a homogeneous polynomial in K with deg{q) < k, then q G 
TL = ©;^gA ^A(<y5), where A„ is the set of eigenvalues ofTn- 

Proof. Let s < k, and let Qs and be as before. Since Lp has accuracy k, and 
s < K, all components of Qb (in fact all components of Qs) are in S{ip), and satisfy 





. 


. 


\ 


1 


(3 . 


. 








. 


• /? 





V 


. 


. 1 





Qnix). (5.6) 



If we denote by (x) the first coordinate of (x) we see that Qg (x) is actually of 
class n{A,P, 1). Hence, by Proposition Ol Qb{x) = r$, where P„y G Ker(r„ - 
/?/). This means that (3 is also an eigenvalue of r„ and the theorem follows. □ 

The following corollary imposes conditions on the eigenvalues of T„, under the 
hypothesis of accuracy. 



Corollary 5.1. Assume that ip has accuracy k and that {ip(- — k)}ki£T ire linearly 
independent. Let Xi,...,Xd be the eigenvalues of A (counted with multiplicity). If 
rj = (yj-, . . . , Y^), then [rj°']\^\^s "^re eigenvalues ofTn, for s = 0, 1, . . . , k — 1. 

Proof. Let Ai, . . . , Ad be the eigenvalues of A. By |CHM98j . [A"]|q,|=s are the eigen- 
values of Also, recall that since A is invertible, A^g] is also invertible and 
(A[s])'^ — A^^[gy So the eigenvalues of ^~^[s] are [r7"]|a|=s. We have already 
proved that if (p has accuracy k and s < k, then every eigenvalue of ^^^[s] is also 
an eigenvalue of r„. So the result follows. □ 

6. 1-DIMENSIONAL EXAMPLES 

We conclude the paper by exhibiting two examples of (A, fc)-homogeneous func- 
tions, associated to scaling functions in dimension 1. The higher dimensional ex- 
amples can be obtained in a similar way. 

6.1. Daubechies D4. Daubechies wavelets, are those refinable functions of iV co- 
efficients, that are orthogonal and provide the highest order of accuracy possible. 
(Note that the splines do not form an orthonormal base). 

D4 is the refinable function that satisfies the refinement equation of 4 coefficients: 



Di{x) = ^^D4{2x) + ^^D4{2x-l) + ^^Di{2x-2) + ^^D4{2x-3). (6.1) 
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Figure 1. Daubcchies D4 with the homogeneous functions. 



Figure 2. ScaUng function for coefficients 1/3, 2/3, 2/3, 1/3 with 
Homogeneous functions - /ii, /i2, /13. /13 is a 2-homogeneous func- 
tion 



D4 has accuracy 2 (it reproduces the constant and the linear functions). 

In this case the matrix T has eigenvalues 1, \ and cq = ^'^^ (see ICHMOSp . So a 
basis for span{D4(a;), r'4(2:-|-l), D4^{x+2)}^fz[o^i\ is also given by spanjl, x, hca{x)}x£[a.i\ 
where hc^ is the homogeneous function associated to cq. 



6.2. (A, 2)-Homogeneous function. In the previous example, we obtained a local 
basis of span{/(a;), /(a;-|- 1), /(x + 2)}^g[Q just by using 1-homogeneous functions. 
The following example is to illustrate, that even in the simple case of only 4 coef- 
ficients, it may be necessary to use homogeneous functions of order bigger than 1 . 
Consider the function: 

f{x) = i/(2x) + ^/(2.x - 1) + '^f{2x - 2) + i/(2x - 3). (6.2) 

It can be shown that / has accuracy 1, and the eigenvalues of T are {1, |}. So in this 
case, span{/(a;),/(a; + 1), f{x + 2)}^g[o,i] = span{l, /i{i/3_i}(a;), /i{i/3_2}(a;)}:jg[o,i]' 
where /i{i/3.i} is a 1-homogeneous function corresponding to the eigenvalue 1/3, 
and /i{i/3,2} is a 2-homogeneous function corresponding to the eigenvalue 1/3 (see 
Figure |2Jl. 
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